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SUMMARY 

An experimental investigation was conducted in a heated-air facility in order to 
determine the effects of a variation in uncooled pipe -inlet diameter and the associated 
contraction area ratio on the heat -transfer distributions in 30° and 60° half -angle of 
convergence nozzles having the same throat diameter and curvature. Tests were con- 
ducted at a nominal stagnation temperature of 970° R (539° K) and nominal stagnation 
pressures of 160 and 300 pounds per square inch (110 and 207 N/cm^) absolute. Experi- 
mental heat -transfer coefficients were compared to values predicted by two interrelated 
boundary -layer calculation techniques and a pipe -flow (Nusselt number) type of correla- 
tion. 

Pronoimced increases in the experimental nozzle peak heat-transfer coefficients 
accompanied a reduction in contraction area ratio, other conditions remaining fixed. 

For the same contraction area ratio and stagnation conditions, a pronotinced increase in 
peak heat -transfer coefficient was coincident with an increase in nozzle convergence 
angle. 

All the prediction techniques failed to yield consistent agreement with the experi- 
mental peak heat -transfer coefficients; that is, in some cases all methods underesti- 
mated the peak values of heat-transfer coefficient, but in other cases an overprediction 
was evident. The overpredicted heat transfer was usually characteristic of the configu- 
rations having the largest contraction area ratio (*«19). 


INTRODUCTION 

The design of regeneratively cooled nozzles for hi^ performance rocket systems 
requires a comprehensive knowledge of the gas-side convective heat transfer. Engineer- 
ing techniques for predicting the heat transfer often consist of using correlations which 
have been tailored to the experimental heat -transfer distributions corresponding to a 
certain class of nozzle geometries and flow environments (ref. 1). These correlations. 


which are usually premised on flat -plate, zero pressure-gradient heat transfer, consider 
primarily the mass -flux variation through the nozzle. The effects of flow acceleration 
and turning are completely neglected. Hence, it is not surprising that a simple corre- 
lation equation failed to predict the heat transfer in a hydro gen -oxygen rocket (ref. 2) 
and heated air nozzles (refs. 3 and 4). 

In reference 5, a Nusselt number correlation was employed in the prediction of heat 
transfer to 30° and 60° half-angle of convergence nozzles operating with heated air. 
Calculated throat heat -transfer coefficients for the 60° half -angle of convergence nozzle 
agreed reasonably well with the experimental results; however, a gross overprediction 
(by nearly a factor of 2) resulted when the correlation was applied to the 30° half -angle 
of convergence nozzle. 

Another recent study, which qualitatively illustrates the inadequacies inherent in a 
correlation approach, is described in reference 6. Results obtained in a center-body 
experiment at sonic conditions revealed a large difference* in friction coefficient for the 
same throat Reynolds number based on diameter. This difference in friction coefficient 
was obtained by altering the nozzle geometry. Through Reynolds analogy, a similar 
effect on the heat transfer would be expected. 

A fundamental method of calculating the nozzle heat -transfer coefficients is pre- 
sented in reference 4. This method involves the simultaneous solution of the integral 
momentum and energy equations in conjunction with Coles friction law and a modified 
von Karman momentum-heat analogy. Application of this boundary -layer theory to air 
nozzle, data in references 3 to 5 yielded better agreement with experimental peak heat- 
transfer coefficients than predictions based on correlation techniques. The radical 
difference in convergence angle of the various nozzles provided a rather severe test for 
the boundary-layer analysis since the relative flow histories in terms of the turning and 
pressure gradient were appreciably different. 

The present experimental study was conducted to further assess the influence of 
inlet conditions on the heat transfer in 30° and 60° half -angle of convergence nozzles; 
however, in contrast to the experiment of reference 5, the diameter rather than the 
length of the uncooled pipe inlets was altered. Inlets having diameters of 6. 50, 3. 07 
and 2. 47 inches (16. 51, 7. 80, and 6. 27 cm) corresponding to contraction area ratios of 
approximately 18. 9, 4.2, and 2. 7, respectively, were used in the investigation. All 
tests were conducted with air at a nominal stagnation temperature of 970° R (539° K). 
Nominal stagnation pressures were 160 and 300 pounds per square inch (110 and 207 
N/cm ) absolute. Experimental nozzle heat -transfer coefficients are compared to 
calculated values based on (1) the boundary -layer analysis of reference 4, (2) an energy 
calculation (ref. 7) which was derived as a limiting case of the theory of reference 4, 
and (3) a Nusselt number correlation. 
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area 

diameter of insulating void surrounding heat -flux meters 
constants 

exponent of Reynolds number in skin-friction law 

skin -friction coefficient 

low-speed adiabatic skin-friction coefficient 

specific heat of air at constant pressure 

local diameter 

film -temperature function 

heat -transfer coefficient 

enthalpy 

enthalpy difference, i^^^ - i^ 

enthalpy difference, Iq " V 

thermal conductivity of Inconel 

thermal conductivity of air 

Mach number 

Nusselt number 

interaction exponent 

pressure 

Prandtl number 

local heat flux 

Reynolds number 

radius 

Stanton number 
distance along wall 
temperature 

temperature on heat -flux meter 
velocity at edge of boimdary layer 
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X axial coordinate measured from nozzle throat 
X length of inlets 

Y distance along heat -flux meter measured from gas-side wall 
z axial coordinate measured from nozzle entrance 
z axial distance from nozzle entrance to throat 
a constant in skin -friction law 

/3 angular position of nozzle instrumentation (tables I and n) 

A temperature thickness of boxmdary layer 
6 velocity thickness of bovmdary layer 
e contraction area ratio 
6 momentum thickness of boimdary layer 
[i viscosity 

p gas density at edge of boxmdary layer 
cp energy thickness of boundary layer 
Subscripts: 

ad adiabatic wall condition 
D based on diameter 
i based on enthalpy 
in incompressible 
p peak value 
r reference condition 

s static condition at edge of boimdary layer 
t throat 
w wall condition 
9 based on momentum thickness 
cp based or energy thickness 
0 stagnation condition 
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APPARATUS 


The nozzle heat -transfer facility comprises the following components: (1) a heat 
exchanger, (2) plenum and bypass flow control system, (3) imcooled pipe inlet, (4) water- 
cooled nozzle, and (5) an exhaust system. The facility, which was the same as that of 
references 3 and 5, is shown schematically in figure 1. 


Heat Exchanger 

The hot gas ejected from a jet -engine burner section was used to heat a series of 
1 -inch -diameter (2. 5-cm-diam) Inconel coils through which passed high-pressure dry 
air. The dry air was heated to a temperature of about 970° R (539° K) at pressures of 

o 

160 and 300 poimds per square inch (110 and 207 N/cm ) absolute. 

Plenum and Bypass Flow Controller 

The plenum consisted of a 6 -foot-long (1. 83-m-long) chamber having entrance and 
exit diameters of 14. 0 and 8. 4 inches (35. 6 and 21. 3 cm), respectively. A Pitot pres- 
sure probe having an inside diameter of 0. 186 inch (0. 472 cm) was positioned on the 
centerline of the plenum. This probe provided a reference pressure for the manometer 
system used to measure nozzle static pressures. 

A flow bypass manifold, coupled to the downstream end of the plenum, was used to 
control the velocity profile upstream of the pipe inlets by removing the boundary layer 
along the plenum wall. Bleed flow rates were adjusted by means of a pneumatically con- 
trolled valve operating in a critical flow mode. 


Uncooled Pipe Inlets 

Uncooled (adiabatic) pipe inlets having inside diameters of 6. 50, 3. 07, and 2. 47 
inches (16. 51, 7. 80, and 6. 27 cm) and lengths of 17. 0, 18. 0, and 18. 1 inches (43. 2, 

45. 7, and 46. 0 cm), respectively, were used in this investigation. These inlets pro- 
vided nozzle contraction area ratios of approximately 18.9, 4.2, and 2. 7, respectively. 
Details of the inlet geometries are shown in figure 2. All three Inlets were used in con- 
jimction with a 30° half -angle of convergence nozzle; however, only the 6. 50- and 
3. 07-inch -diameter (16. 51- and 7. 80-cm-diam) inlets were tested with a 60° half -angle 
of convergence nozzle. All the inlets were machined from AISI 304 stainless-steel pipe 


to a wall thickness of 0. 25 inch (0. 64 cm) and an inside surface finish of 64 micro- 
inches (1.63 Mm) (rms). A thin film of silicone rubber was applied to the downstream 
end of the inlets to minimize the heat conduction from the water-cooled nozzles. Sili- 
cone rubber was also used to seal the 3. 07- and 2. 47 -inch -diameter (7. 80- and 6. 27-cm- 
diam) inlets in the supporting chamber as noted in figures 2(b) and (c). 


Water-Cooled Nozzles 

Two water-cooled conical nozzles, each having a nominal throat diameter and throat 
radius of curvature of 1. 5 inches (3. 8 cm) were used in this investigation. The nozzles 
had a half -angle of convergence of 30° and 60° and a half -angle of divergence of 15° (here- 
after these nozzles will be identified by their half -angles). The nozzles were machined 
from AISI 304 stainless -steel forged billets and had a wall thickness of 0. 5 inch (1. 3 cm). 
Cooling water was directed from the downstream to the upstream end of the nozzles by 
means of Incite shrouds. 

Coordinates for the 30°-15° and 60°-15° nozzles are presented in tables I and n, 
respectively. The expansion ratios for the 30°-15° and 60°-15° nozzles, given by the 
ratio of exit to throat area, were 25. 3 and 3. 32, respectively. The exact contraction 
area ratios corresponding to the various diameter pipe inlets were 18. 98, 4.23, and 
2. 74 for the 30°-15° nozzle and 18. 80 and 4. 19 for the 60°-15° nozzle. 


Exhaust System 

The nozzles were connected to a 24-inch -diameter (61-cm-diam) exhaust duct by a 
thin flexible flange or metallic bellows which yielded to thermal expansion. The exhaust 
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pressure, nominally 2 pounds per square inch (1. 4 N/cm ) absolute, was sufficient to 
prevent nozzle flow separation as confirmed by the measured static pressure at the last 
station in the nozzles. 


INSTRUMENTATION 

Local heat -transfer rates and wall static pressures were measured at stations 2 to 20 
in the 30°-15° nozzle and at stations 2 to 16 in the 60°-15° nozzle as noted in tables I 
and n, respectively. Nozzle wall temperatures were obtained at each of these stations. 
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Pressures 


Nozzle wall static -pressure taps having a diameter of 0. 031 inch (0, 079 cm) were 
coupled to manometers which contained mercury, acetylene tetrabromide, and dibutyl 
phthalate, each selected to give maximum sensitivity in a given Mach number range. In 
the inlet and nozzle entrance region where differences between the static and total pres- 
sure were small, the manometers were connected in a differential fashion and were 
referenced to the total pressure. This manometer coupling technique provided sufficient 
accuracy to detect the presence of an adverse pressure gradient in the entrance region. 
Manometer reference pressures were measured with Bourdon tube gages. Fltdd temper- 
ature corrections were applied in the reduction of the manometer data. 


Heat- Flux Meter 

Steady-state measurements of the gas-side wall temperature and local heat -transfer 
rates were obtained by means of an Inconel plug -type heat -flux meter described in 
reference 8. The diameter of the plug, shown in figure 3, was the same for the two 
nozzles. Three Chromel-Alumel 0.003-inch (0.008-cm) wire thermocouples were spot- 
welded to the 0. 12 5 -inch -diameter (0. 318-cm-diam) Inconel plug in an inert gas environ- 
ment. The thermocouple spacing was determined to the nearest 0.001 inch (0. 003 cm) 
by a microscope. The heat-flux meters were installed with a push fit at the gas side of 
the nozzle and a high temperature O-ring seal at the water-cooled side. A stagnant air 
column surrovmding the plug provided the thermal insulation necessary for one- 
dimensional heat conduction through the shaft. Each heat -flux meter was located 180° 
from the static pressure tap in the same axial plane as noted in tables I and n for the 
30°-15° and 60°-15° nozzles, respectively. 


DATA REDUCTION 

Lxal Heat-Flux and Wall Temperature 

Temperatures measured on the Inconel heat -flux meter were incorporated in the 
Fourier heat conduction equation to determine both the gas-side wall temperature and the 
local heat -transfer rate per unit area. This equation is 


q = -K(t) 

dY 
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where the thermal conductivity K of the Inconel can be represented to within about 
±2 percent over the desired temperature range by the following linear relation: 


K = 1. 307x10"'^ t(°F) + 1. 704x10''^ Btu/(sec)(in. )(°F) 


or 


K = 1. 758x10"^ t(°K) + 8. 243X10“^ W/(cm)(°K) 

Integration of the Fourier equation yields a temperature distribution along the Inconel 
plug of the following form: 

-qY +%t 

The xmknowns q and were determined by the simultaneous solution of two equations 
containing the measured temperature t and the corresponding location Y at two of the 
three stations along the heat -flux meter. Wall temperatures were computed by setting 
Y = 0. 

Measurement of three temperatures along the meter shaft provided three possible 
combinations for computing the nozzle wall temperature. In approximately 95 percent 
of the readings, the average value of the three wall temperatures was used in the subse- 
quent calculation of the local heat-transfer coefficient. Defective thermocouples 
neces;sitated the computation of the wall temperature from a single pair of thermocouples 
in the remaining 5 percent of the readings. 


Error Considerations 

The measured heat fluxes in this investigation, as in reference 3, were estimated 
to be no more than 10 percent higher than the true one -dimensional values. This 
10-percent uncertainty factor, which is due to the void around the heat-flux meter, does 
not apply to the wall temperature. Direct measurements of the peak wall temperature 
by means of Chromel-Alumel thermocouples agreed to within 1 percent of the absolute 
temperature computed from plug measurements. As a matter of consistency, all exper- 
imental heat -transfer coefficients will be presented as computed from the direct 
measurement of one -dimensional heat flux; however, in certain comparisons of the data 
to various predictions, a 10-percent conduction ruicertainty band will be attached to the 
peak heat -transfer coefficient. Experimental inaccuracies in heat flux resulting from 
errors in (a) plug temperatures, (b) thermocouple spacing, and (c) the thermal con- 
ductivity of the Inconel are less than ±3 percent and, therefore, for purposes of this 
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investigation have been considered negligible. Errors in experimental heat flux resulting 
from (1) the relative difference in the thermal conductivity of the stainless -steel nozzle 
and the Inconel plugs and (2) the effects of wall curvature have also been neglected. 


PROCEDURE 

Tests were conducted at nominal stagnation pressures of 160 and 300 poimds per 
square inch (110 and 207 N/cm^) absolute at a stagnation temperature of 970° R (539° K). 
While approaching the desired operating conditions, the bleed flow rate was adjusted 
mitil the plenum boundary layer was removed as confirmed by a mean velocity profile 
obtained with a hot wire upstream of the 6. 50-inch-diameter (16. 51-cm-diam) pipe inlet. 
The same bleed flow rate was established for all inlet configurations. 

All temperature and pressure data were obtained in a period of 15 minutes. Tem- 
peratures were recorded 10 times by an automatic voltage digitizer in order to verify the 
presence of steady-state operating conditions and to permit the averaging of small 
recording errors. The Bourdon tube reference pressure gages were read, and the ma- 
nometers were photographed 2 to 4 times in the 15 -minute recording interval. 


EXPERIMENTAL RESULTS 
Nozzle- Pressure Distribution 

Experimental and one -dimensional static- to total-pressure ratios P^/Pq for the 
nozzles are tabulated in table m. Also shown are the pressure ratios for the imcooled 
pipe inlets; however, static pressures were measured only in the 6. 50-inch-diameter 
(16. 51-cm-diam) inlet. One -dimensional values of pressure ratio were used in calcula- 
ting the boundary -layer development in the 3. 07- and 2. 47-inch -diameter (7. 80- and 
6. 27-cm-diam) inlets. 

The measured pressure ratios deviated considerably from the one -dimensional 
values in the throat region of the nozzle as was also observed in references 9 and 10. 

In addition, static -pressure measurements indicated the presence of an adverse pressure, 
gradient in the corner formed by the intersection of the 6. 50 -inch-diameter (16. 51 cm- 
diam) inlet and the nozzles (table IQ, inlet and station 2). Similar observations were 
made in reference 10 for a cylindrical inlet and conical nozzle intersection having a 
slight radius of curvature rather than a sharp comer as in the present investigation. 

The pressures tabulated in table IQ also indicate that the adverse gradient terminates 
and again becomes favorable between stations 2 and 3 in the nozzle or at a distance along 
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the wall of less than 1.4 inches (3. 6 cm) from the nozzle entrance. In the case of the 
smaller diameter inlets, the presence of an adverse gradient could not be detected with 
the existing instrumentation. 


Heat-Transfer Distributions in the Nozzles 

In this section the experimental heat -transfer distributions, presented in figures 4 
to 8, are discussed. Three distributions, representing predicted results, are also 
shown in this series of figures; however, discussion of the predicted heat transfer will 
be deferred to the succeeding section. The experimental results have been duplicated 
in figures 9 to 13 for purposes of clarity in comparisons with additional theoretical 
curves. 

Experimental heat -transfer distributions will be presented in terms of the heat- 
transfer coefficient based on enthalpy which is given by the following equation: 

^ad " V 


where the adiabatic enthalpy was given by 

*ad = ‘s + - ‘s> 

The Prandtl number was assumed constant at 0. 71. Static and total enthalpies i„ and 

s 

iQ, respectively, were computed using the equations of reference 11 for enthalpy, and 
specific heat in conjimction with an assumed isentropic expansion process. Heat- 
transfer coefficients will be plotted as a function of the axial distance from the nozzle 
entrance z; therefore, caution should be exercised in comparing results for different 
inlet diameters on the basis of the same value of z. The relation between z and the 
axial distance from the nozzle throat x can be obtained from tables IV and V which 
contain the experimental values of h. and T^, respectively, for all the configurations 
and stagnation conditions investigated. Experimental values of h. are plotted as a 
fvinction of axial distance for the 30°-15° nozzle with the 6. 50-, 3. 07-, and 2. 47 -inch- 
diameter (16. 51 -, 7.80-, and 6. 27 -cm -diam) inlets in figures 4, 5, and 6, respectively, 
and for the 60°-15° nozzle with the 6. 50- and 3. 07 -inch -diameter (16. 51- and 7. 80-cm- 
diam) inlets in figures 7 and 8, respectively. Distributions of hj^ are presented for 
each configuration at nominal stagnation pressures of 300 and 160 pounds per square 
inch (207 and 110 N/cm ) absolute corresponding to parts a and b, respectively, in each 
figure. 
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As mentioned previously, the plenum boundary-layer bleed flow rate was the same 
for all inlet configurations. This flow rate was sufficient to remove the plenum boimdary 
layer and establish a uniform velocity profile across the 6. 50-inch-diameter (16. 51-cm- 
diam) inlet. However, in order to establish the effect of possible secondary flows ahead 
of the 3. 07- and 2. 47 -inch -diameter (7. 80- and 6. 27-cm-diam) inlets resulting from the 
method of installation of the smaller inlets (fig. 2), tests were repeated with a zero 
bleed -flow rate. Measured heat -transfer rates with no bleed flow were the same as the 
values obtained with the bleed flow incorporated. Further evidence of this insensitivity 
of nozzle heat transfer to upstream flow conditions is presented in references 3 and 5. 

In all tests a reduction in stagnation pressure for other conditions remaining fixed 

resulted in a lower heat -transfer distribution as expected. If at a given station h. is 

Oft ^ 

assumed to vary as tiie mass flux to a power h.o!(pu) , a decrease in stagnation pres- 
sure of from 300 to 160 pounds per square inch (207 to 110 N/cm^) absolute would 
produce a reduction of 40 percent in h.. Experimentally, the reduction in the peak value 
of h. was between 29 and 47 percent where these limits correspond to the 6. 50-inch- 
diameter (16. 51-cm-diam) inlet coupled to the 30*^-15*^ and 60°-15° nozzles, respec- 
tively, (refer to figs. 4 and 7). 

Comparison of figures 4 to 6 for the 30°-15° nozzle with figures 7 and 8 for the 
60° -15° nozzle reveals that a reduction in the xmcooled pipe inlet diameter produced an 
appreciable increase in the peak value of h.. In the 30°-15° nozzle operating at 
300 pounds per square inch (207 N/cm ) absolute, an increase of 72 percent is the peak 
value of hj was obtained for a reduction in inlet diameter of from 6. 50 to 3. 07 inches 
(16. 51 to 7. 80 cm) (figs. 4 and 5). The maximum value of h. obtained in this study was 
with the smallest inlet which had a diameter of 2. 47 inches (6. 27 cm) (fig. 6). This 
maximum value was 110 percent hi^er than the peak value of hj for the configuration 
having an inlet diameter of 6. 50 inches (16. 51 cm). 

In the 60°-15° nozzle operating at 300 pounds per square inch (207 N/cm^) absolute, 
an increase of 45 percent in the peak value of hj was obtained when the inlet diameter 
was reduced from 6. 50 to 3. 07 inches (16. 51 to 7. 80 cm) (figs. 7 and 8, respectively). 
The same trends of higher heat transfer with lower inlet diameter were apparent for both 
nozzles at the lower stagnation-pressure level. 

Experimental peak values of hj were consistently higher in the 60°-15° nozzle for 
the same inlet and stagnation conditions. A comparison of figures 4 and 7 shows the 
peak value of hj in the 60°-15° nozzle was about 47 percent higher than the value in the 
30°-15° nozzle when operating with the 6. 50 -inch -diameter (16. 51-cm-diam) inlet at a 
stagnation pressure of 300 pounds per square inch (207 N/cm^) absolute. Similar differ- 
ences in the 30°-15° and 60°-15° nozzle peak values of h. were apparent for the 3. 07- 
inch-diameter (7. 80-cm-diam) configuration and also for the lower stagnation-pressure 
operating condition. 
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COMPARISON OF EXPERIMENTAL AND PREDICTED HEAT- 
TRANSFER COEFFICIENTS 


In this section, experimental heat -transfer distributions are compared to the results 
based on (1) a Nusselt number correlation, (2) the boundary -layer theory of reference 4, 
and (3) an energy calculation method generated from the boundary -layer theory. In fig- 
ures 4 to 8, results are given for the Nusselt number correlation as well as the boundary- 
layer theory in which the diabatic skin-friction coefficient was assumed equal to the 
adiabatic value. In figures 9 to 13, the experimental heat -transfer distributions are 
repeated and are compared to predictions based on the energy calculation method and 
the botmdary -layer theory in which the diabatic skin-friction coefficient was assumed 
equal to the product of the low-speed adiabatic value and a film -temperature correction 
factor. Examples of the effects of the static-pressure distribution and interaction expo- 
nent on the heat transfer calculated from the boundary -layer theory are given in fig- 
ures 14 and 15, respectively. Finally, in figure 16, the experimental and the predicted 
peak heat -transfer coefficients are compared on the basis of nozzle contraction ratio. 


Nusselt Number Correlation 

The experimental heat -transfer coefficients are compared to a simple pipe-flow type 
of correlation herein referred to as a Nusselt number correlation given by 

NUj. = 0. 026 Re^- ® Pr^/^ 

where the subscript r indicates that the flow properties are evaluated at the Eckert 
reference enthalpy condition (ref. 12). The Prandtl number Pr was assumed constant, 
0. 71, and the Reynolds number Re^j ^ based on the local diameter of the nozzle is 
given by 


Re 


D,r 




Heat -transfer coefficients were calculated from the Nusselt number according to the 
following equation: 
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hi= 


k^Nu^ 
r r 

C D 


Heat -transfer distributions based on this Nusselt number correlation are presented 
in figures 4 to 8. The maximum deviation from the experimental peak value of hj 
occurred with the 30°-15° nozzle operating with the 6. 50 -inch -diameter (16. 51-cm-diam) 
inlet at a stagnation pressure of 300 pounds per square inch (207 N/cm ) absolute. In 
this case, the Nusselt number correlation overpredicted the peak value of h. by nearly 
84 percent (fig. 4(a)). At the lower stagnation pressure, the Nusselt number correlation 
overpredicted the peak value of h. by about 55 percent (fig. 4(b)). An appreciable over- 
prediction in the peak value of h. can also be noted in figure 7 for the 60°-15° nozzle 
operating with the same inlet. For the other nozzle -inlet combinations, presented in 
figures 5, 6, and 8, the agreement with the experimental peak values of h. was much 
better with a maximum overprediction of 6 percent and a maximum imderprediction of 
15 percent (figs. 5(a) and 8(a), respectively). 


Boundary- Layer Theory 

Experimental values of h-^ have been compared to predictions based on the turbulent 
boimdary -layer theory of reference 4. The method essentially involves the solution of 
the integral momentum and energy equations and employs Coles friction law, 1/7 -power 
profiles for velocity and temperature difference, and the von Karman form of Reynolds 
analogy. Reference 4 should be consulted for details of the assumptions Incorporated 
in this method. Two basic options, provided in this boundary -layer program, concern 
the selection of an interaction exponent and method of evaluating the diabatic skin -friction 
coefficient. 

Inter action exponent option . - The first option concerns the selection of an exponent n 
called the interaction ejqjonent which is used to relate the Stanton number for unequal 
momentum and energy thicknesses to that for equal thicknesses by means of a factor 
icp/Q)^. In the present investigation the interaction exponent was assumed to be zero 
imless specified. 

Ski n -friction option . - The second option deals with the method of evaluating the 
diabatic skin-friction coefficient Cj. One of the methods of determining Cj was to 
assume that it was equal to the adiabatic skin -friction coefficient c^ obtained when 
Tw = The value of Cj was based on free -stream gas properties. 

A second method of determining c^ incorporates a film -temperature correction 
factor which is applied to the low-speed value of c^ (denoted by c^) according to the 
following relation: 
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The value of Cj or was determined by a method of Coles (ref. 13) which 
gives Cj as a fimction of Mach number M and the momentum thickness Re3molds 
number Re^. Details of the calculation procedure are given in reference 4. 

Bovindary -layer in itialization. - The boundary -layer calculation of the nozzle heat 
transfer reqtdres initialization with respect to the momentum thickness 9 and the ratio 
of thermal to velocity layer thickness (A/6 )^_q. Since the inlets were uncooled, only 
the momentum thickness development in the pipes was considered. In all calculations 
involving a zero interaction coefficient, the heat -transfer is independent of the value of 
^z=0’ however, because of the programming technique 0 ^_q was a reqioired input 
variable. Estimates of the value of ^^=0 obtained from the following equation: 

«z=0 " 

where x was the length of the inlet, and the constant, 0. 04, was obtained from the 
boundary -layer measurements of references 2 and 5. 

The boundary -layer calculation was reinitiated at the nozzle entrance using the calcu- 
lated value of and an assumed value of (a/ 6 )^_q. In order to determine the sensi- 

tivity of the assumed value of (A/6 )^_q on the peak value of hj, a range of values of from 
0. 005 to 1. 0 were incorporated in the program. The peak value of h^ was rather sensi- 
tive to the value of (A/6)2 ,_q in the range of 0. 1 to 1. 0; however, values in this range 
were considered too large in view of the adiabatic inlet condition which caused the thermal 
boundary layer to begin at the nozzle entrance. The value of A/6 at this station was 
expected to be quite small, that is, A « 5. For values of (a/ 6 )^_q < 0, 1, the peak value 
of h. rapidly approached an asymptotic value. A lower value of (A/6 )^_q equal to 0. 01 
was used in the calculations because it ensured convergence to the asymptotic peak heat- 
transfer coefficient within a practical computing time. 

Heat -transfer distributions . - The results of the boundary -layer calculation of nozzle 
heat transfer are presented in figures 4 to 8 for assumed values of (A/6 )^_q equal to 
0.01 and 1.0. In this series of calculations, Cj = c^ corresponds to free-stream gas 
properties. Results are presented for the two values of (A/6 )^_q to indicate the sensi- 
tivity of this assumption on the calculated values of h^; however, as stated previously, 
the value of (A/6 )^_q = 0. 01 is expected to better represent the actual flow. 

An examination of figures 4 to 8 indicates that the calculated heat -transfer distribu- 
tion was higher for the lower of the two assumed values of (A/6 )^_q. Also, a comparison 
of figure 4 with figure 6 for the 30°-15° nozzle and figure 7 with figure 8 for the 60°-15° 
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nozzle reveals that the effect of (A/6 )^_q on the peak value of h. was greater for the 
small diameter inlets (low contraction ratio configurations). In general, the best agree- 
ment between theoretical and experimental peak values of hj was obtained with (A/6 )^_q 
equal to 0. 01. An exception, however, can be noted for the 30°-15° nozzle operating 
with the 6. 50 -inch -diameter (16. 51-cm-diam) inlet (fig. 4(a)). In this case, the pre- 
dicted peak value of hj was nearly 35 percent higher than the experimental value. 

The experimental heat -transfer distributions, which have been duplicated in figures 9 
to 13, are compared to predictions based on the second option of the skin -friction coeffi- 
cient. As noted previously, this option assumes Cj = f(T^/Tg)c£. All calculations were 
based on a value of (A/6)^_g equal to 0. 01. 

Incorporation of this second form of the friction law resulted in a downward displace- 
ment of the previously obtained heat -transfer distributions in which c^ = c^ This 
displacement can be noted typically by comparing figure 4 with figure 9 or figure 7 with 
figure 12 for the 30°-15° and 60°-15° nozzles, respectively. This second form of the 
friction law improved the agreement between experimental and theoretical peak values 
of h. in the case of the 30°-15° nozzle with the 6. 50 -inch -diameter (16. 51-cm-diam) 
inlet (figs. 4 and 9) and in the case of the 60°-15° nozzle with the 6. 50-inch-diameter 
(16. 51-cm-diam) inlet (figs. 7(b) and 12(b)). In all the other configurations, the agree- 
ment degenerated. 


Energy Calculation of Nozzle Heat Transfer 


In reference 7, it was shown that for certain assumptions, including one -dimensional 
flow and a zero interaction exponent, a simple closed form solution of the integral energy 
equation can be generated. As shown in the appendix, the assumption of one -dimensional 
flow is unnecessary, and a simplified solution similar to that of reference 7 can be ob- 
tained. The expression for Re^ (eq. (A6)) is 


(Re^Mr Ai) 


1-c _ (1 


- c)q! / ^r Ai 

1.8 / p Ai 

«'z=0 


( cp, r' 

(p 


(Re^ ^lr Ai) 


1-c 


1 + 



dz 


+ (Re^ Mr 


where a and c are constants in the friction law (eq. (A4)). 

A further simplification can be made in the expression for Re if the additional 
assumptions of an isothermal wall, an enthalpy recovery factor of 1. 0, and constant 
viscosity are made. The equation for Re^ (eq. (A7)) then becomes 
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(Re r) 

cp 


1-c _ (1 - c)a I 

1.8 / 

«'z=0 


Pj. , \l-c 

-± S^ii_(Re r)^ ^ 

P (p ^ 


IM- 

\dzy 


0.5 


dz + (Re ^ r)3~n 
' cp 'z=0 


The expression for hj (eq. (A8)) in which a = 0. 0258 and c = 0. 25 is 


h. = 0. 0143 pu Re“°- 

1 rp (pi T 

r 


which can also be expressed in terms of the Nusselt number (eq. (A9)) as follows: 


Nu^ = 0. 0143 Re°-J^ ^ 

This expression for the Nusselt number differs from the value in the simple pipe-flow 
type equation by the factor 0.438 (^j/D)”^’^^. 

The results of the energy calculation of nozzle heat transfer obtained from the 
variable temperature and isothermal wall expressions (eqs. (A6) and (A7), respectively) 
are presented in figures 9 to 13. In these figures, the isothermal wall designation also 
implies a recovery factor of 1. 0. A comparison of the distributions of heat -transfer 
coefficients obtained from the two forms of the simple energy calculation shows that the 
isothermal wall assumption is good for the class of nozzles and the stagnation conditions 
investigated. The isothermal wall assumption generally produced less than 10 percent 
higher peak values of h^ than the variable temperature wall form of the equation as 
noted typically in figures 9 and 12 for the 30°-15° and the 60°-15° nozzles, respectively. 

The results of the energy calculation parallel those of the full boimdary -layer theory. 
Peak values of h^ calculated by the energy method, ranged from 16 to 23 percent higher 
than the values based on the boimdary -layer theory. The similar trends in the results 
are not surprising since the first-order difference between the simplified and the full 
boimdary -layer theory containing a zero interaction exponent are merely the friction law 
and nonunity Prandtl number correction factor. 


Influence of Nozzle-Pressure Ratio on the Calculated Heat Transfer 

In practical applications of the boimdary -layer theory, it is frequently desirable to 
use one -dimensional pressure ratios to describe the mass-flux variation in the nozzle. 

In table IH, pronounced deviations from the one -dimensional pressure ratios can be noted 
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in the form of an adverse gradient in the nozzle entrance and lower values of static pres- 
sure in the throat region. The effect of the pressure distribution on the calculated heat 
transfer thus becomes a primary consideration. 

A comparison of the measured pressure ratio Pg/^o ^ 50-inch-diameter 

(16. 51-cm-diam) inlet and station 2 in the 30°-15° and 60°-15° nozzles (table m) indi- 
cates the presence of an adverse pressure gradient. Since the magnitude of the pressure 
difference appeared to be greatest in the 60°-15° nozzle, this configuration was selected 
for a study of the influence of the adverse pressure gradient on the heat -transfer distribu- 
tion calculated from the boundary -layer theory of reference 4. The calculated heat- 
transfer distribution (fig. 7(a)) is compared in figure 14 to the distribution obtained when 
one -dimensional nozzle -pressure ratios were incorporated in the boundary -layer theory. 
The greatest influence of nozzle -pressure ratio can be noted at the nozzle entrance where 
the adverse gradient was present. Here the calculated value of h^ based on one- 
dimensional pressure ratios was nearly 4 times as great as the value based on the exper- 
imental pressure distribution. It is probable that this large difference in calculated heat 
transfer is partially a reflection of the sensitivity of the mass fliix to the experimental 
pressure ratio. In progressing downstream, the difference in calculated values of h. 
rapidly diminishes. The peak value of hj based on the experimental pressure distribu- 
tion was only about 5. 0 percent higher than the value based on a one -dimensional pres- 
sure distribution. In most practical applications, this difference in the peak value of h. 
is insignificant thus suggesting that in the absence of an experimental or axisymmetric 
pressure distribution, one -dimensional values may be used without appreciably influenc- 
ing the accuracy of the calculated heat transfer in the throat region. 


Variation of the Interaction Exponent 

In reference 4, the limits for the interaction exponent are given as 0 and 0. 25. The 
effect of various interaction exponents on the calculated heat transfer is shown typically 
in figure 15. Heat-transfer distributions, corresponding to interaction coefficients of 0, 
0. 10, and 0. 25, have been calculated for the 30°-15° nozzle operating with the 6, 50-inch- 
diameter (16. 51-cm-diam) inlet at a stagnation pressure of 299. 5 povinds per square inch 

9 

(206, 5 N/cm ) absolute. The calculations were performed using the adiabatic form of 
the friction coefficient corresponding to free-stream gas properties. A pronounced in- 
crease in peak heat -transfer coefficient accompanied the increase in interaction exponent. 
The peak value of h. based on an interaction exponent of 0.25 was 65 percent higher than 
the value obtained with a zero interaction exponent. In this case, the zero interaction 
exponent produced the best agreement between the experiment and theory. In cases 
where the theory vmderpredicted the experimental peak value of h. (figs. 5, 6, and 8) 
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incorporation of a nonzero interaction exponent would effect only a sli^t improvement 
in the agreement; therefore, within the scope of the present tests, the assumption of a 
zero interaction exponent appears reasonable. 


PEAK HEAT-TRANSFER DEPENDENCE ON CONTRACTION AREA RATIO 

In figure 16, experimental and calculated peak heat -transfer coefficients h. p are 
plotted as a fimction of contraction area ratio e^. The four experimental fairings are 
indicative of the heat -transfer variation with contraction area ratio for the 30°-15° 
and 60°-15° nozzles, each operating at a stagnation temperature of about 970° R 
(539° K) and nominal stagnation pressures of 160 and 300 poimds per square inch (110 

and 207 N/cm ) absolute. A 10-percent experimental uncertainty band has been applied 
to the data for, as mentioned previously, heat conduction resulting from the presence 
of an insulating void arovmd the heat meter could effect an increase in the measured heat- 
transfer rates. 

In figure 16, a substantial increase in the experimental values of h. with reduced 

P 

values of e^, (smaller inlet diameter) can be observed for both nozzles and stagnation 

pressure levels. In the 30°-15° nozzle operating at the hi^-pressure level (fig. 16(a)) 

an increase of 72 percent in h. was obtained (assuming the same relative experimen- 

P 

tal errors) for a reduction in of from 19. 0 to 4. 2. The maximum value of h. , 

c i,p’ 

obtained with the minimum value of (2.7), was about 110 percent higher than the value 
of h. „ for e equal to 19.0. A similar increase in h. „ can be noted at the lower 
pressure level (fig. 16(b)). 

Results for the 60°-15° nozzle operating at the high stagnation -pressure level 
(fig. 16(a)) indicate that an increase in h. „ of 45 percent occurred when the value of 
was reduced from 18. 8 to 4. 2. The corresponding increase in h. at the low 

^ 1, P 

stagnation-pressure level was 55 percent (fig. 16(b)). Again, these numbers are prem- 
ised on the assumption of the same relative experimental errors. 

A pronounced difference between the calculated and the experimental variations of 
h. „ with can be noted in figure 16. Values of h. „ calculated from the Nusselt 

X>P c i,p 

number correlation were totally invarient with e_ and convergence angle for given 

stagnation conditions; however, a pronounced increase in h. with decreasing e was 

1 , p c 

obtained experimentally. The constant values of h. calculated from the correlation 

P 

equation are the result of nearly equal throat mass fluxes and diameters for all values of 
e . The influence on the calculated values of h. resulting from differences in the local 
reference temperatures was negligible. In all cases except for the 30°-15° nozzle having 
a contraction area ratio of 4. 2, a pronounced overprediction (by nearly a factor of 2 at 
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II 


the high value of e„) or modest underprediction can be noted for the Nusselt number 
correlation. 

The results of the boundary -layer theory shown in figure 16 indicate the experimen- 
tally observed trend of increasing h. with decreasing e ; however, the rates at 

1, p C 

which h. „ increases, denoted by the slopes of the curves, are not consistent with the 
P 

measurements. Also, for equal values of e_ the bovindary -layer theory fails to predict 
correctly the difference in h. in the two nozzles operating at 300 pounds per square 

9 P 

inch (207 N/cm ) absolute. At the lower stagnation pressure of 160 pounds per square 

p 

inch (110 N/cm ) absolute, the experimental difference in h. for the two nozzles at a 

^3 P 

given value of is within the theoretical predictions only if an uncertainty factor of 
10 percent is applied to the 60°-15° nozzle data. On the basis of absolute rather than 
relative levels of heat transfer, the largest discrepancy between experimental and theo- 
retical values of h. _ occurred in the case of the 30°-15° nozzle having a contraction 

1, p 

area ratio of 19. 0; however, this difference from the experimental value was much less 
than the corresponding difference arising from the Nusselt number correlation. In view 
of the uncertainties associated with the theoretical assumptions, initialization technique, 
and experiment, the boundary -layer theory is considered to show the effect of contraction 
ratio reasonably well. 

The bovindary -layer predictions of h. shown in figure 16 are essentially represen- 

1? P 

tative of the values obtained by the simple energy calculation technique; however, for the 
purpose of clarity, the results from the energy calculation are not presented in this 
figure. An indication of the ability of the energy calculation to comprehend the contrac- 
tion ratio effect on h^ „ can be obtained by comparing, for instance, figure 9 with 

**^3 P no 

figure 11, which correspond to the 30 -15 nozzle having contraction ratios of 19. 0 and 
4.2, respectively. 


SUMMARY OF RESULTS 

An experimental investigation has been performed in order to assess the effect of a 
variation in uncooled pipe inlet diameter and the associated contraction area ratio of the 
heat transfer in conical nozzles having 30° and 60° half-angles of convergence. Addi- 
tionally, three calculation methods were applied in an attempt to predict the heat transfer, 
with particular emphasis on the throat region of the nozzles. These methods included 
(1) an integral boundary -layer analysis, (2) an energy calculation technique derived from 
the boundary-layer theory, and (3) a Nusselt number pipe-flow type of correlation. 
Measurements were obtained at nominal stagnation pressures of 160 and 300 pounds per 
square inch (110 and 207 N/cm^) absolute. All tests were performed at a nominal stagna- 
tion temperature of 970° R (539° K). The principal results of this investigation are as 
follows: 
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(1) Substantial increases in the peak heat -transfer coefficient accompanied a reduc- 
tion in contraction area ratio (smaller inlet diameters) when imcooled pipe inlets are 
used with cooled conical nozzles. 

(2) The simple Nusselt number correlation, in which the heat -transfer coefficient is 
primarily a function of the Reynolds number based on the local diameter, was only able to 
predict the peak heat transfer for certain combinations of inlets and nozzles. In some 
cases, the correlation equation produced nearly a 100 percent overprediction of the peak 
heat-transfer coefficient. The correlation method did not comprehend the experimentally 
observed convergence angle and contraction area ratio effects on the nozzle heat trans- 
fer. 

(3) For the conditions investigated, the integral boundary-layer theory yielded a fair 
prediction of the nozzle heat transfer. A comprehension of convergence angle and con- 
traction area ratio effects was evident although not to a degree consistent with the exper- 
imental results. In general, the integral boimdary -layer approach was better than the 
Nusselt number correlation. 

(4) An approximate formulation of the boundary-layer theory, in which the energy 
equation was solved exclusive of the momentum equation, produced results consistent 
with those obtained with the general boundary -layer theory. 


CONCLUDING REMARKS 

The observed changes in the nozzle throat turbulent heat -transfer rates accompany- 
ing the alterations in inlet diameter or nozzle convergence angle are expected to be 
predominantly due to changes in the acceleration history of the thermal boundary layer 
which commenced at the nozzle entrance. In reference 5, experimental results indicated 
a negligible influence of inlet momentum history on the throat heat flux; however, unlike 
the momentum consideration, the inlet thermal history is expected to contribute to the 
turbulent transport of heat in the nozzle. Therefore, caution should be exercised in 
extrapolating the present results to configurations in which cooled Inlets are employed. 

Since predictions of the nozzle heat transfer from the energy calculation method and 
the integral boundary-layer theory of reference 4 were in good agreement, the simpler 
energy method is recommended for estimates of nozzle heat transfer or thermal 
boundary-layer characteristics. However, it is important to recognize that in applica- 
tions of the energy method, a knowledge of the initial thermal boundary-layer thickness 
is required. In configurations amenable to this type of initialization, the energy method 
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is recommended in preference to the Nusselt number correlation for calculations of the 
nozzle heat transfer. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 22, 1967, 

129-01-09-06-22. 
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APPENDIX - ENERGY CALCULATION OF NOZZLE HEAT TRANSFER 


The assumptions presented in reference 4, when combined with (1) a zero interaction 
exponent, (2) an approximation for the nonunity Prandtl number term in the modified 
von Karm^ form of the Reynolds analogy and (3) a Blasius friction law in which Cj is 
evalviated at ^ = 0 are sufficient to provide a simple closed form solution of the energy 
equation (ref. 7). 

In reference 4, the Stanton number was computed from a modified von Karm^ form 
of Reynolds analogy given by 



(Al) 


where the interaction coefficient (cp/d)^ is a correction factor for unequal momentum and 
energy thicknesses. For the special case of n = 0, the integral energy equation can be 
completely divorced from the momentum equation thus providing the basis for the closed 
form solution of the nozzle heat transfer. The equations for the nozzle heat transfer 
(which differ slightly from those of ref. 7) are derived herein. 

The integral energy equation for axisymmetric flow can be written as follows: 


9__ = i In pur Ai 9? = ^ (A2) 

rpucp Ai ^ ds T(p \Ai/ 

where Ai = iad " V ^ ^ ^0 " V’ 

The denominator of equation (Al) involves a fimction of Pr and c^ which, for a 
nozzle, has a weak axial dependence. This denominator was assumed constant at a mean 
value of 0.9. Incorporation of this constant in equation (A2) yields 


j 

— In pur Ai cp = 
ds 


1.8 cp 



(A3) 


The Blasius form of the skin-friction coefficient, evaluated on the basis of cp = 9, will 
be incorporated in equation (A3). This friction law can be e:q)ressed as follows: 
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where a and c are constants, and Cj is the incompressible value of the friction 
coefficient. As noted in reference 14, the compressible form of the friction law can be 
written in terms of the Eckert reference temperature as 


Cf = 


(‘^f,in) 


where, in this analysis. 


Therefore, 



= a Re 


c 

(p,r 


Cf - 



(A4) 


where a = 0. 0258, and c = -0.25. Substitution of equation (A4) into equation (A3) yields 
the final differential equation for cp, which is 

T /Re \ ^ 

{pwc<p Ai)^"^^ -- In pur<p Ai = — [ £i£) (pur Ai)^"*^ (A5) 

ds Tj. 1.8\ (p / Ai 

Equation (A5) can be readily integrated to give the following equation for the Re dis- 
tribution in the nozzle; 


(Re^ Mr Ai) 


1-C _ (1 - c)0! 

1.8 


f 

I p Ai 

‘'z=0 




(p 


(Re^ Mr Ai) 


1-C 


l + (- 

\dz/ 


10. 5 


dz 


+ (Re^ Mr Ai)J“J ‘ (A6) 

It should be noted that the integrand in equation (A6) is not a function of cp. 

In the special case of an isothermal wall, an assumed enthalpy recovery factor of 
1. 0, and constant viscosity, equation (A6) reduces to the following relation: 
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(Re r)^-^ = 

^ 1.8 



r) 


1-c 


1 + 



dz + (Re 

cp 


r) 


1-c 

z=0 


(A7) 


When equation (A6) or equation (A7) are solved for Re , the local heat -transfer 
coefficient can be determined from the assumed form of Reynolds analogy given by 


h. = 0. 0143 pu — Re"®’ 

1 rp (p,r 


(A8) 


Equation (A8) can be expressed alternately in terms of the Nusselt number as follows: 

v-0.25 


NUj. = 0. 0143 Re®' "J? Pr^. 


& 


(A9) 
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TABLE m. - STATIC- TO TOTAL -PRESSURE RATIOS 


(a) Nozzles 


Sta- 

tion 


30° -15° 

nozzle 

60°-15° nozzle 

Area 

ratio, 

A/Aj. 

Pressure ratio, Pq/Pq 

Area 

ratio, 

A/A^ 

Pressure ratio, P.g/PQ 

Experi- 

ment 

One- 

dimensional 

Experi- 

ment 

One- 

dimensional 

2 

17. 548 

0. 99946 

0.99921 

14.663 

0. 99996 

0.99891 

3 

11.648 

. 99847 

.99827 

9.352 

. 99907 

. 99731 

4 

6.952 

. 99565 

. 99509 

7. 769 

. 99826 

.99610 

5 

5.591 

.99325 

. 99241 

6.299 

.99691 

. 99405 

6 

4.395 

. 98887 

. 98768 

3. 322 

.9811 

. 97825 

7 

3.328 

.9804 

. 97830 

2.373 

.9560 

.95631 

8 

2.410 

.9624 

.95768 

1.341 

.8227 

. 84342 

9 

1.392 

.8399 

. 85703 

1.023 

.5680 

.62815 

10 

1.024 

.5733 

.63117 

1.000 

.4681 

. 52828 

11 

1.000 

.4600 

, 52828 

1.030 

.3383 

.41399 

12 

1.013 

. 3733 

. 44998 

1.090 

.2612 

. 33605 

13 

1.065 

.2785 

. 36165 

1.172 

.2222 

. 27469 

14 

1.156 

.2084 

. 28457 

1.349 

.1765 

. 19898 

15 

1.348 

. 1820 

. 19938 

1.876 

. 1103 

. 10530 

16 

1.873 

.1125 

. 10556 

2.883 

.0545 

. 05050 

17 

3.669 

.0375 

.03417 




18 

8.391 

.0094 

. 00951 




19 

15.070 

.0038 j 

. 00397 




20 

24. 071 

. 0024 

. 00200 i 





(b) Uncooled pipe inlets 


Inlet d] 

lameter 

Nominal 

Pressure ratio, Pg/PQ 

in. 

cm 

area 

ratio, 

A/Aj 

Experiment 

One- 

dimensional 

”~6.To~ 

16.51 

18.89 

0.99932 

0.9993 

3.07 

7.80 

4.21 


.9866 

^2. 47 

6.27 

2.74 


.9683 


denotes use with 30^-15^ nozzle only. 















HEAT -TRANSFER COEFFICIENTS 


Tq, 970® R (-539° K).] 
15® Nozzle 















Inlet diameter, 

3. 07 in. (7.80 cm) 



Inlet diameter, 2.47 in. (6.27 cm) 


[~Axial distance, z 

Stagnation pressure, P 

[) 

Axial distance, z 

Stagnation pressure, Pq 


in. 

cm 

299.7 

206.6 

160.2 

110.4 

in. 

cm 

299.9 

206.8 

157.9 

108.9 



psia 

N/cm^ 

psia 

N/cm^ 



psia 

N/cm^ 

psia 

N/cm^ 



Heat -transfer coefficient based 



Heat -transfer coefficient based I 




on enthalpy, 





on enthalpy, hj 




lb 

g 

lb 

g 



lb 

g 

Ib 

g 



(in. ^)(sec) 

(cm^)(sec) 

(in. ^)(sec) 

(cm^Ksec) 



(in. ^)(sec) 

(cm^)(sec) 

(in. ^(sec) 

(cm^(sec) 

















































0. 215 

0.546 

36. ^10"^ 

25. 5X10‘^ 

21. 3x10'^ 

15.0x10"^ 







.565 

1.435 

49.0 

34.5 

25.8 

18.1 

0.133 

0.338 

43. 7X10’^ 

30. 7X10’2 

25. 5X10"^ 

17.9x10"^ 

1.062 

2.697 

70.0 

49.2 

38.9 

27.3 

.630 

1.600 

71.4 

50.2 

43.8 

30.8 

1.500 

3. 810 

70.8 

49.8 

44.6 

31.4 

1.068 

2.713 

86.5 

60.8 

50.5 

35.5 

1.675 

4. 255 

65.7 

46.2 

38.5 

27.1 

1.243 

3. 157 

72.3 

50.8 

44.8 

31.5 

1.805 

4. 585 

69.5 

48.9 

41.1 

28.9 

1.373 

3.487 

67.6 

47.5 

43.7 

30.7 

1.930 

4.902 

57.2 

40.2 

35.2 

24.7 

1.498 

3.805 

64.0 

45.0 

38.0 

26.7 

2.067 

5.250 

50.1 

35.2 

35.6 

25.0 

1.635 

4.153 

53.8 

37.8 

34.5 

24.3 

2.309 

5.865 

45.2 

31.8 

27.1 
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(a) 6.50-Inch (16. 51 -cm) inside diameter Inlet used with 30° and 60° half-angle of convergence nozzles. 



(b) 3.07-Inch (7.80-cm) inside diameter inlet used with 30° and 60° half-angle of convergence nozzles. 
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(c) 2.47-lnch (6.27-cm) inside diameter inlet used with 30° half-angle of convergence nozzle. 
Figure 2, - Uncooled AISI 304 stain I ess -steel pipe inlets. 







Heat -transfer coefficient based on enthalpy, hj, g/(cm'^)(sec) 



CD 


^ (a) Stagnation temperature Tq, 969.9° R (538.8° K); stagnation 

pressure Pq. 299. 5 pounds per square inch (206. 5 Nlcm^) 

.SJ absolute. 



Axial distance along nozzle, z, in. 



0 4 8 12 16 20 24 28 

Axial distance along nozzle, z, cm 


(bl Stagnation temperature Tq, 970.6° R (539.2° K); stagnation 
pressure Pq, 157.8 pounds per square inch (108.8 N/cm^) 
absolute. 

Figure 4. - Experimental and predicted heat -transfer distributions In 
30°-15° nozzle with 6.50-inch (16.51-cm) inside diameter by 17.0- 
inch (43.2-cm) long uncooled pipe inlet. 
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(a) Stagnation temperature Tq, 965. 3° R (536. 3° K); stagnation pressure 
Pq, 299. 7 pounds per square Inch (206.6 HIcnr) absolute. 



Axial distance along nozzle, z, cm 


(b) Stagnation temperature Tq, 964. 1° R (535.6° K); stagnation pressure Pq, 
160.2 pounds per square inch (110.4 HIcnr) absolute. 


Figure 5. - Experimental and predicted heat-transfer distributions in 30°-15° 
nozzle with 3.07-inch (7.80-cm) inside diameter by 18.0-lnch (45. 7-cm) 
long uncooled pipe inlet. 










Heat -transfer coefficient based on enthalpy, hj, g/(cm^)(sec) 



Axial distance along nozzle, z, cm 

(b) Stagnation temperature Tq, 967.6° R (537.6° K); stagnation pressure Pq, 
157.9 pounds per square inch (108.9 HIcnr) absolute. 

Figure 6. - Experimental and predicted heat-transfer distributions in 30°-15° 
nozzle with 2.47-inch (6.27-cm) inside diameter by 18. 1-inch (46.0-cm) 
long uncooled pipe inlet. 
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0 2 4 6 8 10 12 

Axial distance along nozzle, z, cm 

(b) Stagnation temperature Tq, 975.2° R (541.8° K); stagnation pressure 
Pq, 157.3 pounds per square inch (108.4 N/cm^) absolute. 

Figure 7. - Experimental and predicted heat-transfer distributions in 
60°-15° nozzle with 6. 50-inch (16. 51-cm) inside diameter by 17. 0- 
inch (43.2-cm) long uncooled pipe Inlet. 



GO 

CO 



0 2 4 6 8 10 

Axial distance along nozzle, z, cm 

(b) Stagnation temperature Tq. 968. T R (537. 8° K); 
stagnation pressure Pq, 159.4 pounds per square 
Inch (109.9 N/cm^l absolute. 

Figure 8. - Experimental and predicted heat -transfer 
distributions in 60°-15° nozzle with 3.07-inch 
(7.80-cm) inside diameter by 18.0-inch (45.7-cm) 
long uncooled pipe inlet. 
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0 4 S 12 16 20 24 28 

Axial distance along nozzle, z, cm 

(b) Stagnation temperature Tq, 970.6° R (539.2° K); stagnation 
pressure Pq, 157.8 pounds per square inch (108.8 N/cm^) 
absolute. 

Figure 9. - Experimental and predicted heat-transfer distributions 
in 30°-15° nozzle with 6.50-inch (16.51-cm) inside diameter by 
17.0-inch (43.2-cm) long uncooled pipe inlet. 
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0 4 8 12 16 20 24 28 32 36 

Axial distance along nozzle, z, cm 


lb) Stagnation temperature Tq, 964. 1® R (535.6° K); stagnation pressure Pq, 
160,2 pounds per square inch (110.4 HIcur) absolute. 

Figure 10. - Experimental and predicted heat-transfer distributions in 30°-15° 
nozzle with 3.07-inch (7.80-cm) inside diameter by 18.0-inch (45. 7-cm) long 
uncooled pipe Inlet. 
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(b) Stagnation temperature Tq, 967.6° R (537.6° K); stagnation pressure Pq, 
157.9 pounds per square inch (108.9 N/cm^) absolute. 

Figure 11. - Experimental and predicted heat-transfer distributions in 30°-15° 
nozzle with 2.47-inch (6.27-cm) inside diameter by 18. 1-inch (46.0-cm) 
long uncooled pipe inlet. 



0 2 4 6 8 10 12 

Axial distance along nozzle, z, cm 

(b) Stagnation temperature Tq, 975.2° R (541.8° K); 
stagnation pressure Pq. 157. 3 pounds per square 
inch (108.4 N/cm^) absolute. 

Figure 12. - Experimental and predicted heat -transfer 
distributions in 60°-15° nozzle with 6. 50-inch 
(16.51-cm) inside diameter by 17.0-inch (43.2-cm) 
long uncooled pipe inlet. 



Axial distance along nozzle, z, cm 

(b) Stagnation temperature Tq, 968. 1 ° R (537.8° K); 
stagnation pressure Pq, 159.4 pounds per square 
inch (109.9 N/cm^) absolute. 

Figure 13. - Experimental and predicted heat -transfer 
distributions in 60°-15° nozzle with 3.07-inch 
(7.80-cm) inside diameter by 18. 0-Inch (45.7-cm) 
long uncooled pipe inlet. 
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Figure 14. - Predicted heat-transfer distributions in 
60°-15° nozzle with 6.50-inch (16.51-cm) inside 
diameter by 17.0-lnch (43.2-cm) long uncooled 
pipe inlet using one-dimensional and measured 
static pressures. Stagnation temperature Tq, 
967.9° R (537. 7° K); stagnation pressure Pq, 

300. 2 pounds per square inch (207.0 N/cm^) 
absolute. Theory (ref. 4; c^ = Cj g^). 
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0 4 8 12 16 20 24 28 

Axial distance along nozzle, z. cm 

Figure 15. - Influence of Interaction coefficient on predicted heat-transfer 
coefficients for 30°-15° nozzle operating with 6.50-inch (16. 51-cm) 
inside diameter by 17. 0-inch (43. 2-cm) long uncooled pipe inlet. Stag- 
nation temperature Tq, 969.9° R (538.8° K); stagnation pressure Pq, 
299.5 pounds per square inch (206.5 N/cm^) absolute. 


Nozzle 



Figure 16. - Variation of peak heat -transfer coefficient 
with contraction area ratio. Stagnation temperature 
To,-970°R(-539°K). 
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^'The aeronautical and space activities of the United States shall he 
conducted so as to contribute to the expansion of human knowU 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof/^ 

— ^National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 
importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Scientific and technical information generated 

imder a NASA contract or grant and considered an important contribution to 
existing knowledge. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information derived from or of value to NASA 
activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
Utilization Reports and Notes, and Technology Surveys. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington^ D.C. 20546 


